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ABSTRACT 

The Jeans analysis is often used to infer the total density of a system by relating the 
velocity moments of an observable tracer population to the underlying gravitational 
potential. This technique has recently been applied in the search for Dark Matter in 
objects such as dwarf spheroidal galaxies where the presence of Dark Matter is inferred 
via stellar velocities. A precise account of the density is needed to constrain the ex- 
pected gamma ray flux from DM self-annihilation and to distinguish between cold and 
warm dark matter models. Unfortunately the traditional method of fitting the second 
order Jeans equation to the tracer dispersion suffers from an unbreakable degeneracy 
of solutions due to the unknown velocity anisotropy of the projected system. To tackle 
this degeneracy one can appeal to higher moments of the Jeans equation. By intro- 
ducing an analog to the Binney anisotropy parameter at fourth order, /?' we create 
a framework that encompasses all solutions to the fourth order Jeans equations. The 
condition = /(/3) ensures that the degeneracy is lifted and we interpret the sepa- 
rable augmented density system as the order-independent case (3' = (3. For a generic 
choice of /3' we present the line of sight projection of the fourth moment and how it 
could be incorporated into a joint likelihood analysis of the dispersion and kurtosis. 
The framework is then extended to all orders such that constraints may be placed to 
ensure a physically positive distribution function. Having presented the mathematical 
framework, we then use it to make preliminary analyses of existing data leading to 
interesting results which strongly motivate further study. 
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1 INTRODUCTION 

The favoured ACDM model of cosmology is consistent with 
a large invisible non-baryonic component of matter. To in- 
fer its existence astronomers thus look for the gravitational 
effect of its significant mass upon luminous tracer objects or 

for the observable products of DM annihilation and/or decay 

i it ~\ | 1| — — i J 

such as gamma rays (Gunn et al. 1978; Stcckcr 1978) which 

have been used in recent searches fe.g lAbdo et al.ll2010T ) for 
dark matter. In both instances the density distribution of 
the system is critical with the Earth-incident flux of anni- 
hilation products dependent not only on model -dependent 
prop erties derived from particle physics (see e.g. iPieri et al.l 
2009) but also on the square of the density distribution of 
dark matter within the source. This is encoded in what is 
known as the astrophysical J-factor which can be written 
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|Walker et al.ll201ll ). 
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where d is the distance to the source, Pdm(f) is the local den- 
sity of dark matter and ftnt is the integration angle which 
is related to a given solid angle of the source via Afi = 
27r(l — cos 9i n t)- The quadratic dependence upon the density 
introduces a very significant and DM model independent 
contribution to the flux which makes the choice of a stro- 
physical source critically dWalker fc Penarrubial l201ll) im- 
portant for optimising DM searches. Though one might ex- 
pect that the galactic center would provide the strongest sig- 
nal, the strong and chaotic astrophysical backgrounds make 
it arguably less favourable than dwarf spheroidal galaxies 
(dSphs) of the local group which have b een identified as 
having a large mass-to-light discr epancy (lAaronsonl 1983) 
suita ble for dark matter searches ([Lake! 1 19901 : lEvans et alJ 
2004). This, in conjunction with their relative proximity to 
earth, makes them natural laboratories for DM and in recent 
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years it has been possib le to obtain samples of stellar posi- 
tions and velocities (e.g lWalker et al" 2009a) that are large 
enough for statistical treatment. Since the typical angular 
resolution of gamma ray telescopes is larger th an the angu- 
lar si ze on the sky of dwarf spheroidal galaxies (lAbdo et al.l 
l2010h . it turns out that the J-factor is not extremely sensi- 
tive to the distribution of dark matter in the core of dwarf 
spheroidals, although in the event of a signal being observed 
we would ideally want a better indication of the J-factors 
than the range of current estimates which vary over about 
an order of magnitude. 

Another reason why it is important to study the centre 
of dwarf spheroidals is to probe another aspect of dark mat- 
ter, namely its primordial velocity. In the cold dark matter 
hypothesis the kinetic energy of dark matter at the start of 
structure for mation is some very small fraction of its rest 
mass energy l|Hofmann et alj|200ll ; [G reen et al1l2004h and 
the smallest structures above this free streaming scale are 
the first to form. In mo dels of hot dark matter (see e.g. 
iDoroshkevich et al Jl980h dark matter begins completely rel- 
ativistic and the largest structures form first. A (rather finely 
tuned) compromise b etween these two extremes is the idea 
of warm dark matter l|Bond et al.ll 19821) where dark matter 
is not created with highly relativistic velocities, but never- 
theless with significant velocities, meaning that the normal 
growth pattern of cold dark matter proceeds only above a 
length scale related to the free streaming length correspond- 
ing to the initial velocity. This idea has been invoked to 
explain the lack o f predicted satellites of the Milky Way 
l|Moore et al.ll 19991 ) as well as some interpretations of tracer 
populations in dwarf spheroidals wherein it is argued that 
dark matter halos possess a significan t core, possibly due t o 
some inherent initial kinetic energy (jGilmore et al.l [2007f ) . 
At the same time, the importance of the role of baryons 
upon dark matter density in the core of halos is becoming 
increasingly clear |Governato et al1l2012l) . Whatever the un- 
derlying physics, it is clear that we would like to be able to 
interpret the stellar velocity dispersion in such objects more 
effectively. 

To infer the DM density from the kinematic data the 
Jeans analysis is used to relate the joint distribution of tracer 
stars' positions and velocities to the underlying potential of 
the system. Traditionally the second order Jeans equation 
|Binnev fc Tremainel [20081 ) is used to generate the veloc- 
ity dispersions for a set of input parameters including the 
potential which is then fitted to the dSph data with a like- 
lihood analysis by radially binning the line of sight veloci- 
ties for the variance and assuming Gaussianity. It has long 
been known however (|Deiongheil987l ; lMerrittlll987l ) that this 
analysis may not be used to uniquely specify the potential 
for anisotropic systems for which the variances of the radial 
and tangential velocity components are not equal. As it is 
only possible to observe the projected quantities along the 
line of sight, the intrinsic dispersions of the system are con- 
volved such that there is a degeneracy of indistinguishable 
solutions to the Jeans analysis. Indeed it has been shown 
that the observed line of sight dispersion can be generated 
by any given param eterisation of the anisotropy parameter 
(E vans et ah! 120091 ) thus leaving the potential almost com- 
pletely unconstrained. This is the so-called Jeans degeneracy 
problem which is the main subject of this work. 

A discussion of the higher order moments is presented 



herein with a mathematical description of how they enter 
the Jeans analysis and what assumptions are required to 
ensure that the Jeans degeneracy is solved or at least par- 
tially lifted. This is then placed into the practical context 
of a joint likelihood analysis of the variance and kurtosis in 
dwarf sphero idal galaxies. We evaluate the contribution by 
lLokasI (|2002h in establishing a mod el for the kurtosis that 
may be used to lift the degeneracy l|Lokas et al.ll2005l ) and 
extend the method to general anisotropy as proposed bv lAnl 
(|2011bl ) with the separable augmented density system. 

To simplify the mathematical description of the higher 
order Jeans equations there has been much success in the lit- 
er ature since the ad vent of the augmented d ensity formalism 
bvlDeionghel Il 986). Whilst an application (|Deionghelll987l ; 
iBaes fe van Hesd 120071) of this method has gene rally been 
limited to models l|Plumme J 19151 ; lHernquistlll990l ) with par- 
ticu l arly si mple potential-density pairs, the recent work of 
IXnl (|2011bl ) describes for generic density and anisotropy, a 
separable system that solves the Jeans degeneracy problem 
completely by specifying moments at all orders with the po- 
tential and anisotropy parameter alone. This is however by 
no means a general solution and without a strong physical 
motivation its practical use is limited. The layout of our 
paper is as follows. In section [5] we review the mathemat- 
ics behind the Jeans equation and line of sight calculations. 
To utilise the constraining power of the fourth order statis- 
tics we are thus motivated to provide the full analytical set 
of fourth order solutions which is achieved with the intro- 
duction of an analog to the Binney anisotropy parameter 
at fourth order. This is outlined in section [3] wherein we 
show how to construct a generic model for the projected 
fourth order moment and how one could incorporate it into 
a joint likelihood analysis, extending t o full genera l ity th e 
over-constrained method employed by lLokas et al.l (2005). 
In section [4] we use data from the Fornax galaxy to con- 
struct probability distributions for both the J-factor and the 
density of dark matter at lOOpc from the galactic centre as- 
suming an NFW profile. We contrast directly the results of 
the traditional and joint analysis and by direct comparison 
of fourth order models are able to evaluate the robustness 
of observables to the choice of ft' . Finally we will make some 
concluding remarks and outline our future research program. 



2 PRELIMINARY 



2.1 Moments of the Distribution Function 

In the study of stellar systems, a 6 -dimensiona l function / 
is used to specify the distribution (|jeanslll915l ) of stars in 
position and velocity space. For a spherically symmetric sys- 
tem this is related to the underlying gravitational potential 
$(r) by t he time-independent an d collisionless Boltzmann 
equation (|Merrifield fc Kentlll990l ), 
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dt r dr \ r dr I dv r 

+ -(vacate -vrvg)^- (2) 
r dvg 

1 df 

- -UUIV + V^Vg COt 6)— 

r ovj, 
= 0. 

Multiplying ([2]) by vlv^v"^ and then integrating over all ve- 
locities restates it in terms of its true velocity moments 



vvpvfvf = / «f vf vf f{r, v)d 3 v. 



(3) 



where v(r), as an effective zeroth moment that marginalises 
the distribution function in velocity space, is the local den- 
sity of stars. Due to the spherical symmetry of the system 
it is trivial to show by averaging over the azimuthal angles 
that the odd moments vanish and that many of the true 
even moments are related by constant prefactors. To make 
the notation more compa ct we again follow the example of 
iMerrifield fc Kent! |l990) and introduce the intrinsic mo- 



ments of the tangential velocity vt = (vg + v 



2N1/2 



v^vfvf = ~B(j + -,k + -)w* 



„2(i+fc) 



2''" 1 2'~ T ^ 
where B(x,y) is the Beta function. In the subsequent anal- 
ysis we find that it is not necessary to explicitly refer to the 
true moments at any stage and to simplify the mathematics 
the tangential moments will be used exclusively from here 



2.2 Jeans Equations 

To isolate the second order moments i.e the radial and 
tangential dispersions which in practice have the small- 
est statistical errors, the Boltzmann equation is tradition- 
ally multiplied by v r an d integrated over all velocities 
jBinnev fc TremaindlioOsT ) to give, 

d(ua 2 ) , 2/3 _ 2 , _ d* _ n 



2/3 2 

H vo r + v— — = 0. 

dr r dr 



where v(r) is the local stellar density, $(r) is the gravita- 
tional potential that depends on the total density of the 
system p(r) = v(r) + p dm (r) via, 

AttG 



$( r ) = 



r p(r)dr 



(6) 



and t he Binney anisotropy parameter (jBinnev fc Tremainei 
EH) /3(r), 



2*Ur)' (7) 
measures the deviation of the dispersions from the isotropic 
system (o> = a'g = \a\ ) wherein all directions in velocity 
space are equally probabl^B- F° r a dSph, wh ere the mass- 
luminosity ratios are often greater than 10 (|Mated 1 19981 ) 
the dark matter component is very significant. 



1 We choose for mathematical convenience to adopt the 2 D tan- 
gentia l dispersion rather than the ID employed by e.g lLokasI 
(2002) which accounts for the additional factor of 2. 



To illustrate the higher order analysis we consider first 
the fourth order where multiplying equation ([2]) by v% and 
v r Vg respectively relates the three intrinsic moments at 
fourth ord er v*,vf and v 2 vf by th e two fourth order Jeans 
equations |Merrifield fc Kenj[l99o1 '). 



d(i>vf 
dr 



—^-o 2 — r „ 2 d$ 
vivi H — vv^ + 6va r —— 
r dr 
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(8) 
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(9) 



dtvviv?) 1 — 

dr r r dr 

The a dvent of the augmented density system by iDeionghei 
(|l986h has greatly enhanced the mathematical description 
of the Jeans analysis and it is within this framework that 
the co mplete set of Jeans equations has been presented (|Anl 
l2011ah . 



d(yv r 
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dr 



2 2q+2 



(10) 



The number of equations at 2nth order is therefore n, the 
number of permutations of (p, q) for which p + q = n and 
l<p<n, < g < n. Each of the n moments at 2nth order 
enter the derivative of a corresponding equation apart from 



2.3 Projected Moments 

Unfortunately due to the distant nature of astronomical ob- 
jects the stellar positions and velocities may only be ob- 
served along the line-of sight and rather than the true mo- 
ments and the local density one must instead use the pro- 
jected moments and the surface density profile E as functions 
of the projected radius R. The surface density profile is the 
projection along the line of sight of v{r), 

v[r)r 



E(R) 



R Vr 



R 2 



zdr. 



(11) 



which may be inverted directly via the Abel Inversion for the 
local density. The component of the velocity along the line 
of sight, which for convenience is chosen as the z-direction, 
may be expressed as 



R 2 R 

V\os = v r cos(a) — vg sin(a) =v r \ll— — ^ — v$ — - (12) 

where sin a = R/r is the unobservable depth angle that de- 
termines the extent to which each stellar velocity is radial 
or tange ntial. Finding the projecte d moment at 2nth order 
(see also lDeionghe fc Merrittlll992l ) is akin to averaging over 
the 2nth power of equation (|12p . 



_jAx)r 



£<"(#) = 2/ K cos a -we sin a) 

Jr Vr 2 — R 2 



dr 



2 / VCW.feW 

IR fc-n 



2( „- fc) u(r)r 

Vr 2 - R 2 



and one thus requires all of the intrinsic moments to calcu- 
late the projection with the coefficients, 



R 2 



(14) 



© 0000 RAS, MNRAS 000, 000-000 



4 T. Richardson & M. Fairbairn 



To simplify the projected dispersion the anisotropy parame- 
ter is again introduced in place of the tangential dispersion, 



3 R 2 



R? 



-dr 



(15) 



such that specifying the potential and anisotropy parameter 
is sufficient to first calculate the radial dispersion via equa- 
tion ([5]l and then the projected moment for comparison with 
observation. An explicit calculation of the projected fourth 
moment 



Evf^R) = 2 / Ga.owjS + Ca.iigu? + C 2 ,2v' t 



-.dr 



v(r)r 
Vr 2 - R 2 ' 
(16) 



C 2 .k — 



R \2 

— ~pr) 
(r 2 - R 2 ) 

3 R 4 



k = 
k = 1 
k = 2 



(17) 



recovers the result in iMerrifield fe Kent! l| 199ft ). 



2.4 The Degeneracy Problem 

The normal Jeans analysis, which approximates the distribu- 
tion function by its second order moments, has traditionally 
been the only viable means of modeling the limited sam- 
ple sizes afforded by kinematic surveys of dwarf spheroidal 
galaxies. Un fortunately it has been demonstrated (e.g 
lMerrittlfl987h that the integral equation (|15[) can be highly 
degenerate with no way of distingui shing the entangled 
intrin sic dispersions (although see ijWalker fe Penarrubial 
120111 ) for an interesting recent attempt to break this de- 
generacy). In a typical statistical treatment, where a model 
for the dark matter density is fitted with a hand-picked 
form of the anisotropy, there are numerate parameter sets 
p — {/3(r), &{r)} that yield identical line of sight dispersions 
within statistical errors. As the anisotropy parameter is a 
completely unknown degree of freedom it is impossible to 
uniquely specify the potential with the line-of-sight disper- 
sion alone and such a treatment is liable not only to impre- 
cision but also to inaccuracy. With a recent improvement 
in the data, there has been a renewed interest in the higher 
order moments that may be used to distinguish those param- 
eter sets degenerate at second order. Unfortunately whilst 
the fourth moment is practically within reach , a theoretical 
problem persists. As suggested bv lAnl l)2011al ) we note that 
at each successive order the Jeans analysis introduces n + 1 
moments and only n constraining Jeans equations such that 
the intrinsic moments are not specified by the second order 
parameters p. A minimal requirement to define the system 
at fourth order is to specify one of the fourth order moments 
or, in analogy with Binney's anisotropy parameter, to spec- 
ify the ratio of two. To lift the degeneracy one additionally 
desires that the projected fourth moment depend only upon 
the second order parameters such that a new net constraint 
is added to the system without expanding the parameter 
space. The anisotropy parameter in particular, whilst inher- 
ent to the dispersions, has no direct bearing on the higher 
order moments without artificial insertion. It can be con- 
cluded therefore that to utilise the higher order moments 



one must present a new constraint to the system via a sim- 
plifying assumption or empirical observation that optimally 
ties the higher order moments to the anisotropy parameter 
and thus may be used to constrain it. To lift the degener- 
acy completely on e requires all of the projec ted moments 
which as proved bv lDeionghe fe Merrittl (|l992T ) is equivalent 
to knowledge of the distribution function. 



3 ANALYTICAL SOLUTIONS TO THE JEANS 
DEGENERACY 

In this section we discuss analytical models that place ad- 
ditional constraints on the intrinsic fourth order moments 
for a practical use in tackling the mass-anisotropy degener- 
acy. One pro minent example in the literature is the work of 
iLokad l|2002j ) who assumes a form (|Henonlll973r ) of the dis- 
tribution function f(E, L) = fo(E)L~ 2f3 with a constant ve- 
locity anisotropy that is separable when expressed in terms 
of the specific binding energy E — — $ — \ (y 2 + v 2 ) and the 
specific angular momentum \L\ = rv t . With this assumption 
then without specifying the form of the energy component 
fo(E) it is possible via equation ((3} to calculate the ratio of 
higher order moments with the anisotropy parameter. The 
projected fourth moment is thus shown to be dependent only 
on the second order parameters and is applied to the de - 
generacy problem for the Dra co galaxy (Lokas et al.ll2005l ). 
It has been argued however |Anl 12011a ) that this model, 
whilst relatively clear to interpret, over-constrains the prob- 
lem and may be extended to general anisotropy with the 
separable augmented density model. As such we outline this 
model and provide a formula for the projected moments. 
This system however only provides one solution and whilst 
mathematically elegant it does not yet have a strong phys- 
ical motivation. We therefore provide a framework without 
appealing to the augmented density formalism directly, that 
encompasses all solutions to the fourth order by introducing 
an analog to the anisotropy parameter at fourth order. The 
aim of this framework is both to evaluate the particularly 
convenient separable augmented solution and to facilitate 
more physically motivated models perhaps inspired by em- 
pirical observation. The framework is then extended to all 
orders. 



3.1 The Separable Augmented Density Model 

A full account of the augmented density formalism is not 
necessary for the wider context of this p aper but we di- 
rect the interested reader to the literature l|Deiongheili986l ; 
iDeionghe fc Merritdll992l ; lAnll20Tlal lbl) . Here we give a brief 
outline and then present the main findings relevant to a 
study of dwarf spheroidal galaxies in the simplest possible 
terms. B y careful con sideration of equ ation Q and the Jeans 
theorem (|Jeanslll915l ) it was noted bv lDeionghd (1986) that 
the degeneracy of distribution functions which fit the lo- 
cal density v(r) can be represented by augmenting the den- 
sity into an infinite set of degenerate bi-variate functions 
of the augmented configuration space (^,r) and that pro- 
viding a specific form of the augmented density is equiva- 
lent to providing the distribution function. Here ip is the 
positive binding specific potential defined by ^p- = — 4Jr- 
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Moments are also augmented and to retrieve the observ- 
ables the augmented quantities are deaugmented by tak- 
ing the limit — > *l/(r). The power of the method lies in 
the ease with which one may relate the augmented quanti- 
ties and in augmented configuration space the Jeans equa- 
tions are replaced by the sim ple partial differential equations 
l|Deionghe fc Merritdll992h , 



vv 2p v 2q {^, r) = (2p — l)vu 



2p 



-2 2q 
Vt 



(*,»■) 



(18) 



[r 2q+2 vv 2 "v 2q (*,r)] =(P-|) 



r 2q uv 2 /- 2 v 2q+2 {^,r) 



and the augmented total derivative I Anl (|2011al ). 



dr dr dr 



(19) 



(20) 



To determine the moment ratios at fixed order it is necessary 
only to consider the order preserving equation (|19p such that 
for a system with a separable augmented density, 

v(^,r) = P(^)R{r) (21) 

the anisotropy is dependent only on the radial comp onent 
R(r) and this has been demonstrated by lAnl (|2011al ) who 
shows that 

dlog R 



(22) 



dlog r 2 

An even more profound property of the separable system 
(equation I21[) is that the ratio of moments at a given or- 
der scales only with the constant prefactor inherent to the 
isotropic system. In practice this has the nota ble effect o f 
collapsing all of the Jeans equations to just one l|Anll2011ah . 



d(vv* n ) 2/3 — . „, 

' + — vv 2n + (2n - l)uvr 
r 



d<5> 



, =0 (23) 
dr r dr 

that enables a unique calculation of the radial intrinsic mo- 
ment and thus absorbs the additional degree of freedom in 
the Jeans analysis. In calculation of the other moments one 
may use (lAnlbOllaT ). 



2(n — k) 2k 



Oik 



(24) 



where (o)& = ]~Ii=i( a + i — 1) is the Pochhammer symbol 
that describes the rising sequential product, oto = 1 and 

aq+1 = (q + l-l3)a g + ~^Z (25) 
2 dr 



can be used to iteratively generate all moment ratios noting 
again that one requires only the anisotropy parameter. The 
projected moments then follow trivially, 



= 2/ J2 C - 



'(n-*+§)* 



R 2 



dr (26) 



and we thus demonstrate that the fourth projected moment 
depends only on the second order parameters. In practice 
one must first solve n differential equations to calculate in 
turn each radial intrinsic moment via the Jeans equation 
(|23[) and then n iterations of (|25[) for a theoretical predic- 
tion of the projected moment. As an example consider the 
projected fourth moment for which we require, 



ax = 1-0, a 2 = (l-/3)(2- / 8)-I^. 



(27) 



and the fourth order Jeans equation 

dw£ 2/3 — 2 d<f> 

— ; 1 vv^ + iva r — — 0. 

dr r dr 



Evaluating (I26|l with n = 2 yields after some algebra, 



zdr 



Evf os (R) = 2 g{P,r,R) 



Vr 2 - R 2 



2 p^ + ^ l - p)Ri 



R 4 d/3 
4r 3 dr 



(28) 



(29) 



(30) 



g{f),r,R) . _ 2 

which generalises the result of lLokasI <|2002h and may be 
used to tackle the Jeans degeneracy for an arbitrary radial 
depende nce of the ani s otrop y parameter via the method em- 
ployed in lLokas et al.l (|2005l ). In practice then one must first 
find the radial dispersion via the second order Jeans equa- 
tion |[BJ|, solve equation (]28[) for the fourth radial moment 
and then with the recursive relations (|27|l it is possible to 
evaluate the projected moment with equation ([16}. 

In summary the separable augmented density system 
completely lifts the Jeans degeneracy by providing an in- 
finite set of projected moment equations that one can, in 
principle, calculate from the potential and the anisotropy 
at second order. The origin of the system is however purely 
mathematical and whilst particularly convenient to use it is 
not unique in exhibiting this behaviour. Indeed one could 
artificially introduce any arbitrary relationship between the 
anisotropy parameter and the higher order moments to the 
same effect. This model, which generates only one DF for 
a given set of anisotropy and density parameters, is thus a 
severe restriction to apply to the system without physical 
motivation. The aim then is to provide a framework that 
explores the full range of solutions to the Jeans degeneracy 
problem without restriction from mathematical considera- 
tions. 



3.2 An Extended Model of Anisotropy 

In this subsection we will demonstrate, without appealing 
explicitly to the augmented density formalism, that a unique 
set of of all deprojected moments is specified by introducing 
an analog to the Binney anisotropy parameter at each suc- 
cessive order. Indeed one can represent all possible distribu- 
tion functions in this way and we will show that a distribu- 
tion function can be defined as an infinite set of anisotropy 
parameter analogs {/3n}, 1 < n < oo. If each of these is 
assumed to have a known relation to the anisotropy param- 
eter then the Jeans degeneracy is completely lifted. For a 
particular definition of anisotropy parameters, inspired by 
the separable augmented density system, we provide formu- 
lae for the full set of intrinsic moments and their subsequent 
projected moments. We are then able to interpret the sep- 
arable augmented density system as a particular subset of 
this system for which the anisotropy parameter analogs are 
independent of order {/3 n } = P Vn. 



3.2.1 Fourth Order- 
To measure the anisotropy at fourth order we choose for 
simplicity the adjacent moment ratio to the radial fourth 
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moment and introduce, 



3 v 2 v? 

1 - 



(31) 



which is an analog of the Binney anisotropy parameter 
for the fourth order that measures the deviation from the 
isotropic system (|A8I) where v$ = Substituting this 

parameter into the Jeans equation ||S). 



1 M£ + 3^(7,.— 

dr r dr 



0. 



(32) 



which then, given a functional form for /3' , enables the radial 
fourth moment to be determined uniquely. The choice of 
equation (|3ip is now clear as we note that the separable 
augmented density system thus corresponds to ft = /3'. The 
mixed moment then follows trivially from the definition of 
and the tangential moment may be calculated via equation 
([9]) as follows. Firstly equation (|31[) indicates that, 



d{w^v 2 ) 
dr 



dr dr 



(33) 



which may be simplified further with equation (|32[) by sub- 
stituting in the derivative of the radial fourth moment, 



d(vv?v'j) 2 _ 
dr 3 l P> 



W - , 2 d<S> 

i/v? — 6vo~ r —— 

r dr 



2 dp' - 

3 -dr V V t 



(34) 

Plugging this into the Jeans equation ([9| and rearranging 
for the fourth tangential moment yields, 



-r 4 



r dP' 



= - ( (1 - p'){2 - 0) - - ^ ] vf, + 2(/3' P)ral- 



>d$ 



dr 



(35) 



which shows that the system is uniquely determined by the 
second order parameters plus the anisotropy at fourth order. 
One also notes that in the limit /3' = fi, the tangential mo- 
ment recovers the result vf = ^ct2V$ for the separable aug- 
mented density system. With all the moments one can then 
compute the projected moment from equation (|16[1 which 
after some algebra is, 



Xvf os (R) = 2 



i>(r)r 
Vr 2 - R 2 
(36) 

where g(/3') is adapted from equation (|30|l with /?—>/?'. An 
arbitrary choice of f3' in addition to the second order param- 
eters represents the complete set of fourth order systems. 

Therefore to lift the Jeans degeneracy one must addi- 
tionally impose that /3' = /(/3) where a specification of the 
function f(/3) defines the model. This is the fundamental 
concept of this paper. In section [4] we explore the effect of 
this choice on the key physical observables of the dark mat- 
ter problem. If one can use empirical evidence to find a cor- 
relation between the anisotropy at both orders or indeed a 
well motivated physical argument for the specific form that 
such a relationship would take, then the Jeans degeneracy 
problem can be at least partially lifted with available data. 
To lift the degeneracy completely, all projected moments 
are required and in the following section the framework pre- 
sented above is extended. 



-.dr 



3.2.2 Anisotropy of the Entire System 

To generalise the method in the previous section we again 
choose to define the anisotropy parameter analog at 2nth 
order via the adjacent moment ratio to the radial moment, 



Pn(r) = 1 - (n - ~)^==3- 



(37) 



which measures the departure from the isotropic system, 

— („ _ 1 /r^„,2n-2„2 



= (n - l/2)v 



and we note that Binney's 



anisotropy parameter j3 — j3\ and its fourth order analog 
/3' = P2 are naturally incorporated into the analysis. Sub- 
stituting this into the first Jeans equation at 2nth order, 



d{vv 2n ) 2/3. 



d<b 



dr 



^-vv 2n + (2n - lW™-^— = 

dr 



(38) 



which uniquely determines the radial moment if one assumes 
knowledge of all moments at preceding order. The definition 
of the separable augmented density system is then extended 
to {/3 n } = /3, Vn and thus represents the system for which 
the anisotropy parameters are independent of order. We will 
show that calculation of the other moments then follows via 
a recursive passage through the remaining Jeans equations. 
To generalise the method used to obtain the tangential mo- 
ment in equation (|35[1 we consider the generic Jeans equation 
(1101) and introduce, where convenient, the compact notation, 



2p 2q 

VV T V t — m 



(39) 

and the adjacent moment ratios at order 2n = 2(p + q), 
In m v . a 



f P q + "({M,r) 



2 m p +i, g _i 



(40) 



where we note that by definition /" = 1 — /3 n . The prefactor 
p + 1/2 is included to absorb the order dependence inherent 
to the isotropic system (|A8[) such that we can isolate the 
order dependence of the systems anisotropy. Rearranging 
equation (|38|) for the unknown moment, 



2p-2 2q+2 



r d(yv r ' v 
2 



2pr "t q ) 



dr 



+ (q+ l)uv 2 r P V 2 t q 



I In 2p-2 2qd§ 

+ XP~Tjrvv/ v t q — 



(41) 



such that upon substituting the prior moment ratio (|40[) into 
the derivative, 



dm v 



dr 



P + 



dfq . rndm p+ l }q -i 

mp+1 ' q ~ 1 ~dr~ tq 



dr 



1 dft 2 (fn . 
J^17 + -r {fq ~ q \ 



(42) 



-2f q m p , q -i — 



we obtain the recursive relation for q > 1, 

fg+l 



1 m p _i, g+ i 
v ^ 2' m 



+ 



1 + P q + 
(n-q- 



p.q 

r 1 dK 



l.m„_i,„ d$ 
-) — - — —r — 

2 trip q dr 



(43) 
(44) 



f n 
J'i 



Starting with /" = 1 — j3 n this relation may then be used to 
iteratively calculate the complete set of intrinsic moments 



© 0000 RAS, MNRAS 000, 000-000 



Analytical Solutions to the Mass-Anisotropy Degeneracy 7 



at 2nth order. The first iteration yields 

fS = 2-j3n 



d/3 n 3 TO„_ g _i l8 

- — -r- 



2(1 — p n ) dr 2 m n _ 9 , q dr 1 — /3„ 

(45) 

which we may use to check that with fl = 1 — ft' and 
fi=l — /3, the tangential fourth moment in equation (|35[) 
is recovered via vf = ^fif%Vr- To calculate the projected 
moments we first note that 



q TT 9 f " 

rn n -k,k _ TT m n-j,j _ llj=iJj 

m ",o fJi m n -j + i,]-i {p+\) q 



(46) 



and that in analogy to equation (|25p in the separable aug- 
mented density system we may define, 



*q — I I Jk 

fc=l 



which satisfies the recursive relation 



a q + 1 = Ct q f 



q+1- 



(47) 



(48) 



To check that the system converges to the separable aug- 
mented density in the limit {Pn} = P where the order de- 
pendence is removed, 



fq+l({Pn}) -> fq+l(P) =1+/? + 



2f q dr 



(49) 



we substitute equation (|49p into equation (]48|1 with a q — 
f q a q -i to eliminate f q such that after some algebra, 



= 1 + 



r a q 
2~a 



1 da„ 



da a 



1 da, 



8-1 



a q -i 



2 a 



9-1 



fir 



dr 



Q, + 



9-1 

2~dr~ 



rir 



(50) 



(51) 



which one can prove inductively is satisfied by 

, , . r da„ 
Q9+1 = (<l + k)a q + 2~^T- 

Applying the boundary condition ai = 1 — f3 and choosing 
ao = 1 yields k = 1 — f3 and therefore recovers equation 
Q25p . To generalise the projected moments in equation Q26p 
we simply promote the ratios to a£. Whilst seemingly a 
subtle change, this makes the subsequent calculation con- 
siderably more cumbersome. Though the generalisation at 
fourth order is straightforward the higher order projected 
moments become rather inconvenient for practical use when 
one deviates from the separable augmented density system. 
The inclusion of the higher order analysis is however useful 
to ensure that the set {fin} yields positive moments at all 
orders and thus a physical distribution function. This con- 
straint on the system may be represented simply as, 



fq > 0, V n, q, 



(52) 



which imposes physical constraints on the analogs including 
P n < 1 for q = 1 and for example the sufficient but not 
strictly necessary condition for constant anisotropy /3 n > 
Pn-i corresponding to q = 2. 

3.3 Summary 

By introducing an analog of the Binney anisotropy param- 
eter at fourth order it is possible to represent the complete 



analytic set of proj ected f ourth moments. Not only is the 
result presented in lLokasI |2002i ) adapted in equation (|36l) 
to an arbitrary specification of the anisotropy parameter 
but as a subset of the separable augmented density system 
it is interpreted as the particular case where anisotropy is 
order- independent . To employ the method in lLokas et al.l 
(120051 1 to lift the Jeans degeneracy one may choose any 
arbitrary specification /?' = f(fi) to construct a physical 
model thus removing the constraints upon the distribution 
function imposed by the separable augmented density. Ad- 
ditionally the higher order projected moments are presented 
suc h that c onstraints ()52p analogous to those presented in 
lAnl (|2011al l for the separable augmented density system, 
a q > 0, may be used to ensure that any given model has a 
physical DF. Though the maths underlying the higher order 
moments quickly becomes impractical as one deviates from 
the separable augmented density system, the fourth order 
generalisation remains tractable. With a generic framework 
we facilitate a model with a stronger physical motivation 
and state, in general terms, a condition to lift the degen- 
eracy namely a correlation between the anisotropy param- 
eter analogs. Without a strong physical argument to pro- 
vide this relation analytically we turn to empirical evidence. 
The additional freedom afforded by the anisotropy parame- 
ter analogs provides the means to test models presented in 
the literature. 



4 A PRACTICAL APPLICATION TO THE 
DWARF SPHEROIDAL GALAXY FORNAX 

With the extended anisotropy framework a natural progres- 
sion is to devise a model for f3' with a stronger physical basis 
than the separable augmented density system. Due to the 
completely unknown nature of the anisotropy this is how- 
ever a formidable task and at time of writing no such can- 
didate has been found. To implement the above formalism 
then, we look to dSph data to provide information about 
higher order anisot ropy and as a me ans of eval uating the 
models proposed bv lLokasI |2002i ) and[Xn| (|2011bh . Recently 
the kinematic surveys of dSphs in the Local Group have ex- 
panded the number of velocity measurements to the order 
of 10 3 which reduces the fourth order statistical errors to 
tolerable levels. 



4.1 Data 

In our analysis we use data from the Magellan survey 
l|Walker et alj|2009al ) and in particular the heliocentric line- 
of-sight velocities and RA-DEC coordinates for the Fornax 
galaxy. To remove potential interlopers we use the spectral 
metallicity data provided in 1 Walker et al.l (|2009bh and re- 
move all stars with a membership probability of 95 er cent 
or less retaining a sample of 992 stars. Firstly the Fornax 
positional data were converted to a projected radius from 
the centre of the galaxy assuming a distance of d — 138kpc. 
This was then used to fit the the local density v(r) to a 
Plummer profile by matching the surface profile E(_R) inte- 
grated by the area element dA = 2-uRdR to the number of 
stars at radius R with the result r v — 600Pc. The velocity 
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data was then sorted into 12 radial bins of equal size such 
that there are 12 variance and kurtosis Q measurements. 



4.2 Likelihood Analysis 

Though less favoured than the Einasto jEinasto fc Haudl 
and the more general Zhao density profile (|zhadll996[ ) 
which generalises the slope parameters of the galaxies cen- 
tral cusp, for simplicity at this early stage of our work we em- 
ploy the two paramet er Navarro- Frenk- White (NFW) profile 
ijNavarro et al.lll996h . 



Po 



■z- (l+ — 



(53) 



to parameterise the density of dark matter which is suf- 
ficient to illustrate the statement of concept presented in 
this paper. To model the anisotropy parameter we use 
jBaes fc van Hesell2007r i, 



P(r) = (fin -fa)- 



+ r' 



+ Po 



(54) 



which generically describes anisotropy with a quadratic 
transition about rp that asymptotes at Po for r — and 
Poo for r — Y oo. In addition to the NFW density parame- 
ters {po,r n } the content of the full parameter set depends 
on the nature of the functional form of the anisotropy con- 
sidered at both orders. We consider five permutations (i) 
the traditional second order analysis with a single con- 
stant anisotropy parameter such that the parameters to 
be varied are p — {/3, po,r n } , (ii) a second order analy- 
sis with a varying anisotropy given by equation (|54l) p = 
{Po, Poo, rp, Po, r n}, (hi) a fourth order analysis with param- 
eters equivalent to case (i) for which both a nisotropy pa- 
rameters are constant and equal (|Lokasil2002r ) (iv) a fourth 
order analysis with parameters equivalent to (ii) where both 
anisotropy parameters have the generalised profile but are 
equal, i.e the separable augmented density and (v) a fourth 
order analysis for which both anisotropy parameters are con- 
stant but not necessarily equal p = {/3, P' , po, r n }. 

To fit the anisotropy and density parameters p to the 
empirical data d we employ a Monte-Carlo Markov Chain 
(MCMC) to generate the posterior distribution function 
from the likelihood, 



\nC{d\p) = \ \J2 



+E 



(ftlos Mos) 



(55) 



where N is the number of radial bins, X — {cr los , ki os } are 
the theoretical predictions generated by the parameter set 
p for those quantities evaluated at bin radius Ri, X is the 
empirical measurement of X for the ith bin of data set d 
and ax is the error associated with X. We approximate the 
unbiased sampling errors by, 



N 



2 2 

-<y , 



24iV(iV-2)(iV-3) 
(JV + l) 2 (iV + 3)(jV + 5) 



(56) 



2 The kurtosis is the normalised fourth moment and we note that 
the kurtosis of the normal distribution is 3 



log 10 [l-A] : 


[-1,1] 


loglO r H ■ 


[0,4] 


logio PO ■ 


[-5,3] 


logio : 


[1,5] 



which for a given bin size are generated by considering the 
number of times one must sample from a Gaussian distribu- 
tion to accurately determine its moments. To implement the 
MCMC the Metropolis-Hastings ^Metropolis et al.lll953h al- 
gorithm is used to iteratively sample the parameters or sim- 
ple functions of the parameters for which a uniform sampling 
is desired. As an example consider the parameter P for which 
we choose instead — log 10 [l— P] to sample uniformly over the 
range [ -1,11 i n line with the choi ce for constant anisotropy 
used in lCharbonnier et al.l l|201ll ). This choice equivalent to 
— 9 < P < 0.9 ensures that all ratios of the dispersions are 
uniformly sampled rather than the non-linear phase space 
of p. With this in mind the following priors and ranges are 
adopted, 



(57) 



where A = {P, /?', p , Poo}. To build the MCMC chain a set 
of parameters p' is drawn randomly via a Gaussian distri- 
bution centered at p. The set p' is added to the chain if 
£(d\p') /£(d\p) > U where U is generated randomly from 
the uniform distribution over the range [0,1]. The width of 
the Gaussian known as the temperature of each parameter is 
then adjusted to give an acceptance ratio of approximately 
0.2. 



4.3 Results 

In analysing the fourth order statistics our aims are twofold. 
Firstly we intend to demonstrate the constraining power of 
the fourth order method and by direct comparison with the 
traditional analysis we also look for inconsisten cies. Though 
this h as been achieved in part by the work of lLokas et al.l 
(2005), our second aim, to check the robustness of the fourth 
order model, is new work. Fig. [T] demonstrates the effect of 
introducing the fourth moment to the astrophysical J-factor 
and the DM density profile at lOOpc with Lokas's model. 
Though the width of the distribution is notably reduced for 
the J-factor the improvement upon the precision of the den- 
sity is very small. A noticeable discrepancy is observed be- 
tween the peaks of the J-factor distributions which is relaxed 
when the second order is extended to general anisotropy. 
Whilst one might be tempted to hail this as a success of 
the accuracy introduced by the fourth order statistics, the 
arbitrary nature with which the model P = P' is chosen 
suggests caution and it is entirely feasible that a different 
choice could align the two distributions perfectly. Be that 
as it may, by not fitting to the fourth order at all an as- 
sumption is made to the fourth order without reference to 
the data. In Fig. [2] we demonstrate that if Lokas's model is 
assumed to be true then the second order fit with equiva- 
lently constant anisotropy provides an inaccurate account of 
the kurtosis which decreases the joint likelihood by a factor 
of almost 2. A key question then is whether the uncertainty 
associated with the choice of model at fourth order is sig- 
nificantly less than that which may be alleviated with its 
inclusion. To check the robustness of the Lokas model we 
performed a likelihood analysis for two additional models 
with a simple extension to general anisotropy within the 
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Figure 1. Posterior distributions for a point-like approximation 
to the astrophysical J-factor JTJ with an integration angle of 1° 
(above) and the density of dark matter at lOOpc (below) derived 
from parameters generated via the MCMC for a traditional sec- 
ond order analysis with constant anisotropy (solid), a traditional 
analysis with generalised anisotropy ( I"i4[ ) (dashed) and a joint 
analysis of dispersions and kurtosis with constant anisotropy (dot- 
ted). 
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Figure 3. Distributions as specified in Fig. [T] for joint fits to 
the dispersion and kurtosis with fourth order models /3'(r) = f} 
(dotted), /9'(r) = j3(r) (dashed) and for which /3' is a constant 
that is allowed to vary freely (solid). 
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Figure 4. Best fits to the joint likelihood analysis for the three 
models outlined in Fig. \3\ represented with a consistent colour 
scheme. 



Figure 2. Dispersion and Kurtosis profiles for the three models 
outlined in Fig. [l] The dotted lines show the moments generated 
by the parameters that maximise the joint likelihood function H55H 
for 12 equally sized radial bins of the Fornax data with errors de- 
fined by equation H56I I. For the second order models (solid and 
dashed) the dispersion is generated by parameters that maximise 
the likelihood for a fit to the dispersions only. The resulting kur- 
tosis profiles are produced with the separable augmented density 
model of the fourth moment. 

separable augmented density framework and a model with 
constant anisotropy for which we allowed a constant choice 
to /3' to vary freely. The result of this analysis is shown in 
Figs. [3] and [4] with the surprising result that the additional 
freedom afforded to the system has little impact on the pre- 
diction for the observable quantities shown, indeed drasti- 
cally less than the variation shown from the second order fit. 



Whether this is a property of the data, which being strongly 
platykurtic excludes large portions of the parameter space, 
or the model is unclear and warrants further investigation. 



5 DISCUSSION 

The Jeans analysis is extremely useful in identifying the 
gravitational potential from the line of sight velocities of 
stars moving in that potential and is used to learn more 
about many different kinds of astrophysical objects. The 
information obtained in this way is limited by degenera- 
cies which exist as a result of our ignorance of the velocity 
anisotropy within the stellar tracer populations. By look- 
ing at not only the width (2nd moment of velocity) of the 
stellar velocities but also their kurtosis (ratio between 2nd 
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and 4th moment) we are able to break some of these de- 
generacies to a greater or lesser extent. It has been shown 
that the additional information provided by the higher or- 
der moments can demonstrably tighten the constraints on 
observa bles such as the a strophysical J-factor and the core 
density (jLokas et al.l2005l ). To date this has only been shown 
when an arbitrary and quite constraining choice has made 
to relate the two projected moments to each other. 



In this paper we have presented a new mathematical 
framework for calculating the higher order moments of the 
Jeans equation based upon introducing an analogue of the 
Binney Anisotropy parameter at each higher order and we 
have demonstrated that this completely determines the so- 
lution to the Jeans equation at each order. By studying For- 
nax we demonstrate that a model with constant and order 
independent anisotropy parameters j3' — /3 predicts both 
the J-factor and the density of dark matter at lOOpc from 
the galactic centre with a smaller uncertainty than when 
only the second order moments are taken into a ccount. This 
is ess entially a reproduction of the analysis of (lLokas et al.l 
120051 ) and reaches the same conclusion as reached by those 
authors. We then go on to consider two more additional cases 
which to our knowledge have not been tested against data 
before. The first is the employment of the separable aug- 
mented density, which corresponds to the situation where 
the second order velocity anisotropy /3(r) can be an arbi- 
trary function of radius but where it is equal to its fourth 
order sibling, i.e. P'(r) = /3(r). The second new case we 
consider is where we allow separate anisotropy parameters 
for the second and fourth order moments, although for the 
time being we only consider the case where they are both 
(different) constants. It appears that the predictions for the 
J factor and the core density in the three different cases 
where we have made different assumptions about the sec- 
ond and fourth order velocity anisotropy parameters do not 
vary a huge amount whichever particular set of assumptions 
is adopted - simply by including the higher moments of the 
data, the Jeans analysis seems to become more powerful. 
Whilst a detailed investigation is required to confirm this 
finding it provides strong motivation for further study into 
the higher order Jeans analysis and demonstrates that even 
a fourth order model without a strong physical basis could 
prove more constraining than a second order analysis that 
makes an assumption without reference to the fourth order 
data. 



Ideally we would like to test how constraining the ad- 
ditional fourth order information is in the case where we 
allow the relationship between f3 and j3' to vary to a greater 
or lesser extent. It would also be interesting to try and 
use the results of N-body simulations to motivate physical 
choices for the relationship between the two anisotropy pa- 
rameters. Furthermore, it remains to consider different func- 
tional forms for the dark matter halo such as the Einasto 
profile which can describe both cored and cuspy halos in 
order to see what the higher order information is telling us 
about the shape of the density profile in the centre of dwarf 
Spheroidals. We are working on all these issues and hope to 
present new results in the near future. 
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APPENDIX A: ORDER DEPENDENCE OF 
MOMENT RATIOS IN THE ISOTROPIC 
SYSTEM 

Let's assume that the center of the galaxy lies in the positive 
z-direction for a sphere with corresponding polar angles (0 < 
ei < 7r, < £2 < 2-7r) centered at the stellar position such 
that the radial velocity is treated as v z and the angular 
velocities as v x and v y . We may thus write, 



Vg = V Sin £i COS £2 



Vtf, — v sm ei sm 62 



v r = v cos ei 



(Al) 



(A2) 



(A3) 



where v is the total stellar velocity v 2 — v 2 + v 2 + and we 
note that Vt = v sin t\ . The moments may thus be expressed 
as. 



v 2p v 2q = «2(p+9) CO s 2 P ei sin 2 « ei 



(A4) 



such that to calculate the moment ratios at 2nth order it is 
sufficient in the isotropic case to consider only the angular 
contribution fl — cos 2p ei sin 29 ei . Performing the average 
over all solid angles we integrate dQ = — d(cos ei)de2 with 
the uniform probability distribution P(cosei, £2) = 1/47T, 

Q, = I I P(cos ei, £2) cos 2p ei sin 2q t\ sin e\de\de2 



1 

4tt 
1 



cos 2p ei sin 29 tid cos eicfe 



(A5) 



X 2p (l- X 2 ) q dX. 



The ratio of tangential moments at fixed order is thus, 



m a ,b a 



r(a+|; 



(A6) 



m c , d d!r(c+|) 

where we have used the compact notation (|39[) and the fact 
that for moments of equivalent order a + b — c + d. The 
expression is simplified further for adjacent (a = c — I, b = 
d + 1) moment ratios wherein the shift property of gamma 
functions may be exploited, 

m ™ = -?«U (A7) 
m p +i, 9 _i 2p + I 

Radial moment ratios may then be expressed as the product 
of adjacent moment ratios, 

q-l 

11 _ m "-q+^-' (as) 



n 



L(n- 



h + k) (n - 



2)1 



(A9) 
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and we note crucially that the order dependence of the 
isotropic system is the same as that of the separable aug- 
mented density system. This factor will be used to gener- 
alise the factor of 1/2 present in the second order Binney 
anisotropy parameter. 
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